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ABSTRACT 

The analysis of the 7r7 transition form factor provides severe constraints on the 
pion’s wave function. This information is used to examine charmonium decays into 
two pions critically. It will be argued that the standard perturbative QCD analysis 
of these reactions fails, i. e. the need for additional contributions can convincingly 
be demonstrated. Colour-octet admixtures to the charmonium states are proposed 
as a possible dynamical mechanisms to solve the puzzle. Consequences of the 7r7 
analysis for the electromagnetic form factor of the pion are also discussed. 


1. Introduction 

At large momentum transfer the hard scattering approach (HSA) provides a 
scheme to calculate exclusive processes. Observables are described as convolutions of 
hadronic wave functions which embody soft non-perturbative physics, and hard scat¬ 
tering amplitudes Th to be calculated from perturbative QCD. In most cases only 
the contribution from the lowest-order pQCD approach in the collinear approxima¬ 
tion using valence Fock states only (termed the standard HSA) has been worked out. 
Applications of the standard HSA to space-like exclusive reactions, as for instance 
the magnetic form factor of the nucleon, the pion form factor or Compton scatter¬ 
ing off protons revealed that the results are only in fair agreement with experiment 
if hadronic wave functions are used that are strongly concentrated in the end-point 
regions where one of the quark momentum fractions, x, tends to zero. As has been 
pointed out by several authors (e.g. |@, |^), the results obtained from such wave func¬ 
tions are dominated by contributions from the end-point regions where perturbative 
QCD cannot readily be applied. Hence, despite the agreement with experiment, the 
predictions of the standard HSA are theoretically inconsistent for such wave func¬ 
tions. It should also be stressed that the large momentum transfer behaviour of the 
helicity-flip controlled Pauli form factor of the proton remains unexplained within the 
standard HSA. 

Applications of the HSA to time-like exclusive processes fail in most cases (e.g. 
Gmi -Ptt, 77 —^ VV)- The predictions for the integrated 77 ^ tttt cross-section 
(Icos^l < 0.6) are in fair agreement with the data whereas the predictions for the 
angular distribution fails. Exclusive charmonium decays constitute another class of 
time-like reactions. If the end-point region concentrated wave functions are employed 
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again, the standard HSA provides results in fair agreement with the data in many 
cases. It should be noted that in most calculations of exclusive charmonium de¬ 
cays [Q Ug values of the order of 0.2 — 0.3 are employed. Such values do not match 
with as evaluated at the charm quark mass, the characteristic scale for these decays 
{aairric = 1.5 GeV) = 0.37 in one-loop approximation with Aqcd = 200 MeV). Since 
high powers of a* are involved in charmonium decays a large factor of uncertainty is 
hidden in the predictions. 

Constraining the pion wave function P, |] from the recent precise data on the Try 
transition form factor [0, one observes an order-of-magnitude discrepancy between 
data and HSA predictions for charmonium decays into two pions. In contributions 
from the ccg Fock state are suggested as the solution of this puzzle. I am going to 
discuss these topics in my talk. Also I shall discuss the large momentum transfer 
behaviour of the pion form factor in the light of the new information on the pion’s 
wave function. 

2. The TT-y transition form factor 

The apparent success of the end-point concentrated wave functions, in spite of 
the theoretical inconsistencies, prevented progress in understanding hard exclusive 
reactions for some time. Recently, with the advent of the CLEO data on the Try 
transition form factor 0 , the situation has changed. The leading twist result for 
that form factorf including o^-corrections, reads 



( 1 ) 


/tt is the usual pion decay constant (130.7 MeV) and /ir represents the renormaliza¬ 
tion scale. The function has been calculated by Braaten in the MS scheme. 
(x“^) is the 1/x moment of the pion distribution amplitude, 0, which represents 
the light-cone wave function of the pion integrated over transverse quark momenta, 
k_L, up to a factorization scale, /ij?, of order Q. The distribution amplitude can be 
expanded upon Gegenbauer polynomials, the eigenfunctions of the evolution 

kernel for mesons 



CyH2x - 1 ) ( 2 ) 


^ 4>as[x) 1+ y] B„(^o) I —y- 

n=2,4,... 


where the asymptotic distribution amplitude is (j)As{x) = 6x{l — x). The 1/x moment 
of the distribution amplitude reads 



^ The pion mass as well as the light current quark masses are neglected throughout. 
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Figure 1: The scaled 7 r 7 transition form factor vs. Q^. The solid (dashed) line rep¬ 
resents the results obtained with the modihed HSA using the asymptotic (Chernyak- 
Zhitnitsky) wave function. The evolution of the Chernyak-Zhitnitsky wave function 
is taken into account. The dotted line represents the limiting behaviour Data 

are taken from 0 . 0 - 


The process-independent expansion coefficients Bn embody the soft physics; they 
are not calculable at present, p-o is a typical hadronic scale, actually /iq = 0.5 GeV. 
Since the anomalous dimensions, 7 „, are positive fractional numbers increasing with n 
(e.g. 72 = 50/81) any distribution amplitude evolves into the asymptotic distribution 
amplitude for InQ^ —> cx); higher order terms are gradually suppressed. Hence, the 
limiting behaviour of the transition form factor is 


7T7 


V2U/Q^ 


( 4 ) 


which is a parameter-free QCD prediction |^. As comparison with the CLEO data 
0 reveals, the limiting behaviour is approached from below. At 8 GeV^ the data only 
deviate by about 15 % from (|[) (see Fig.|). In order to give a quantitative estimate of 
the allowed deviations from the asymptotic distribution amplitude one may assume 
that B 2 is the only non-zero expansion coefficient in (^. The truncated series suffices 
to parametrize small deviations. Moreover, it has the advantage of interpolating 
smoothly between the asymptotic distribution amplitude and the frequently used 
Ghernyak-Zhitnitsky distribution amplitude |]^ {B 2 = 2 / 3 ; C 2 ^‘^{^) = 3 / 2 ( 5 ^^ ~ !))• 
For large momentum transfer the assumption on the expansion coefficients is justihed 
by the properties of the anomalous dimensions 7 ^. 

In it is shown that the leading twist, lowest order pQGD result (P nicely hts 
the GLEO data for = —0.39 ±0.05. Using Braaten’s result for 0 that, 

choosing = Q, reads 


10 1 -59/72H2(Q2) 

y 1 ± H2(Q2) 


( 5 ) 
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one finds = —0.17 ± 0.05 from a fit to the CLEO data. Braaten’s analy¬ 

sis is however incomplete in so far as only the corrections to the hard scattering 
amplitude have been considered but the corresponding corrections to the kernel of 
the evolution equation for the pion’s distribution amplitude were ignored. As has 
been shown by Muller recently next-to-leading order evolution provides logarith¬ 
mic modifications in the end-point regions for any distribution amplitude, i.e. for the 
asymptotic one too. An estimate however reveals that the modifications of the evo¬ 
lution behaviour in next-lo-leading order are very small and can safely be neglected. 

To summarize the iA.y form factor requires a distribution amplitude in a leading 
twist analysis that is narrower than the asymptotic one in the momentum transfer 
region of a few GeV^. The Chernyak-Zhitnitsky distribution amplitude is in clear 
conflict with the data and should, therefore, be discarded? 

Recently a modified HSA has been proposed by Botts, Li and Sterman 
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which transverse degrees of freedom as well as Sudakov suppressions are taken into 
account. This approach has the advantage of strongly suppressed end-point regions. 
Hence, the perturbative contributions can be calculated self-consistently in the sense 
that the bulk of the perturbative contribution is accumulated in regions of reasonably 
small values of the strong coupling constant. It is to be stressed that the effects of the 
transverse degrees of freedom taken into account in the modihed HSA represent soft 
contributions of higher-twist type. Still, modihed HSA calculations are restricted to 
the dominant (valence) Fock state. Another advantage of the modihed HSA is that 
the renormalization scale can be chosen in such a way that large logs from higher 
order perturbation theory are eliminated. Such a choice of the renormalization scale 
are accompanied by Ug singularities in the end-point regions which are, however, 
compensated by the Sudakov factor in the modihed HSA. Singularities produced by 
the evolution of the wave function are also cancelled by the Sudakov factor. 

Adapting the modihed HSA to the case of vry transitions, one can write the 
corresponding form factor as [^, |] 

= y dx —-b, ^f) Th (x, b, Q) exp [-S (x, 6, Q)] (6) 


up to 0{as, k^/Q"^) corrections, b is the quark-antiquark separation and is canon¬ 
ically conjugated to the usual transverse momentum k^. The use of the transverse 
conhguration space is mandatory because the Sudakov exponent S is only known 
in that space [0. The Sudakov exponent comprises those gluonic radiative correc¬ 
tions not taken into account in the evolution of the wave function. Th is the Fourier 


^ In |r^ a modification of the pion wave function is proposed where the distribution amplitude 
is multiplied by the exponential exp [—to^o^/x( 1 — x)]. The parameter represents a constituent 
quark mass of, say, 330 MeV. Since the exponential substantially deviates from unity only in the 
end-point regions it leads to a strong additional suppression in the case of the Chernyak-Zhitnitsky 
distribution amplitude ((x“^) changes from a value of 5 to 3.71 at the scale /io). For narrow 
distribution amplitudes {B 2 < 0), on the other hand, the exponential has only a minor bearing on 
the results for 
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transform of the lowest order momentum space hard scattering amplitude. It reads 

Th {x, b, Q) = -^^^0 (Vl -xQbj (7) 

where Kq is the modihed Bessel function of order zero. Due to the properties of the 
Sudakov exponent any contribution is damped asymptotically, i.e. for ln(Q^//rQ) —> 
oo, except those from conhgurations with small quark-antiquark separations and, as 
can be shown, the limiting behaviour (^) emerges, b plays the role of an infrared cut¬ 
off; it sets up the interface between non-perturbative soft gluon contributions - still 
contained in the hadronic wave function - and perturbative soft gluon contributions 
accounted for by the Sudakov factor. Hence, the factorization scale /xf is to be taken 
as 1/6. 

Finally, d'jr is the Fourier transform of the momentum space (light-cone) wave 
function of the pion for which a Gaussian k^-dependence is employed 


/vr f k‘^ ' 

{x,kr,fiF) = ■^^(j)n{x,HF)Nexp ( 


( 8 ) 


Here N = 167 r^a^/(x(l —x)) and, for a distribution amplitude with = 0 for n > 4, 
= l/( 7 r/^Y^ 8 (l -F B 2 ). The 7 r° —>■ 77 constraint []^ is automatically satished for 
that choice of the transverse size parameter Ott. represents a soft wave function, 
i.e. a full wave function with its pertubative tail removed from it. For B 2 = 0 the 
wave function leads to a valence Fock state probability of 0.25 and a r.m.s. radius 
of 0.42 fm. Using the wave function in a modihed HSA calculation, one hnds 
excellent agreement with the CLEO 0] and CELLO data above ~ 1 GeV^ for 
-B 2 (/^o) = —0.006 ±0.014 P, 1^ (see Fig.^). Hence, the asymptotic wave function, i.e. 
the asymptotic distribution amplitude combined with the Gaussian k^-dependence, 
works very well if the modihed HSA is used. 

A similar analysis of the 77 and the ^'7 transition form factors has been carried 
through in 1^. It is important thereby to take into account mass corrections and the 
7]—7]' mixing. The results of that analysis are in excellent agreement with the available 
data including the recent CLEO data [0 . The values of the rj — tj' mixing angle and 
the decay constants are calculated in p to be 9p = —18° ± 2°, = 175 ± lOMeV 

and frj' = 95 ± 6 MeV, respectively. The 7]c'y transition form factor can be analyzed 
in the same manner. Estimates of that form factor can be found in IT? . 


With the advent of the CLEO data the Try transition form factor attracted much 
interest and, besides [|, |[, many papers have been devoted to its analysis elucidating 
various aspects of it [^, |^, 00.01. Particularly interesting is the generalization 
of (|^ to the case of two virtual photons. In the standard HSA and again with = 0 
for n > 4, the Try* transition form factor reads H 


£ ^ 

= a/2 ^ {[1-a;^±2a;ln(n] 


x[l ± 


1 ± 28a; ± 




B2{i^f)] + 10a;lna;i?2(/rF)} 


(9) 
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where uj = Q'‘^/Q‘^. The larger one of the two photon virtnalities is denoted by Q^, 
the smaller one by The factorization scale may be chosen as cts- 

corrections to can be found in [^] and an estimate of power corrections in . 
The treatment of vry* transitions within the the modihed HSA is straightforward 
generalization of (|^) |IB| . 

Interestingly, the form factor still behaves as Q~^ at large Q^. This is to be 
contrasted with the behaviour of the vector meson dominance model |23[ . 

In the limes a; —> 1 @ simplihes to 


F-nry* 


1 + -(1 — a;)(l — 12i?2(/^F)) 


( 10 ) 


3 Q2 

The limiting behaviour of the form factor for a; = 1 which is strictly independent on 
the form of the distribution amplitude, has also been derived from QCD sum rules 
In the triangle diagram is analyzed with the most general form of the irqq 


vertex. The result obtained for F,,^* in that paper is similar to (^) provided B 2 is put 
to zero in (|). The differences between the two results are strongest at cn = 0 (about 
9%) while both the results coincide at a; = 1. 


3. Pionic decays of charmonium 

In view of the results for F^^^ a fresh analysis of the decays Xcj is in order. 

Using the information on the vr wave function obtained from the analysis of F.,^^, one 
hnds the following values for the partial widths 


r(Xco( 2 ) ^ vr+TT-) = 0.872 (0.011) keV (11) 

within the standard HSA |§]. As usual the renormalization and the factorization 
scales are identihed in that calculation and put equal to the c-quark mass. The 
parameter describing the Xcj state is the derivative R'p{0) of the non-relativistic 
cc wave function at the origin (in coordinate space) appropriate for the dominant 
Fock state of the Xcj, a cc pair in a colour-singlet state with quantum numbers 
^s+iLj = ^Pj. rric = 1-5 GeV and, of course, the leading order standard HSA value 
-0.39 for i? 2 (/xo) are chosen as well as R'p{0) = 0.22 GeV^'^^ which is consistent with a 
global £t of charmonium parameters as well as with results for charmonium radii 
from potential models . 

In P the modihed HSA is also used to calculate the Xcj tttt decay widths. 
Taking B 2 = 0 and the other parameters as quoted above, one hnds 


r(Xco(2) ^ 71 + 71 -) = 8.22 (0.41) keV. (12) 

For comparison the experimental data as quoted in and reported in a recent paper 
of the BES collaboration are 

F(xco tt^tt”) = 105 ±30 keV (PDG), 

62.3 ± 17.3 keV (BES), 
r(Xc 2 ^ vr+TT-) = 3.8 ±2.0 keV (PDG), 

3.04 ± 0.73 keV (BES). (13) 
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One notes that both the theoretical results, (ph]) and (|^, fail by at least an order 
of magnitude. To assess the uncertainties of the theoretical results one may vary the 
parameters, rUc, B 2 and Aqcd- However, even if the parameters are pushed to their 
extreme values the predicted rates are well below data. Thus, one has to conclude 
that calculations based on the assumption that the Xcj is a pure cc state, are not 
sufficient to explain the observed rates. The necessary corrections would have to be 
larger than the leading terms. A new mechanism is therefore called for. 

Recently, the importance of higher Fock states in understanding the production 
and the inclusive decays of charmonium has been pointed out |^. It is therefore 
tempting to assume the inclusion of contributions from the |cc 8 (^S'i) 5 f) Fock state to 
exclusive Xcj decays as the solution to the failure of the HSA. The usual higher Fock 
state suppression by powers of 1/Q^ [0 where Q = me in the present case, does 
not appear as a simple dimensional argument reveals: the colour-singlet and octet 
contributions to the decay amplitude behave as 


M 


(c) 


fjf-rn-"’ 


(14) 


The singlet decay constant, represents the derivative of a two-particle coordinate 
space wave function at the origin. Hence it is of dimension GeV^. The octet decay 
constant, /} , as a three-particle coordinate space wave function at the origin, is also 
of dimension GeV^. Since M^f' is of dimension GeV, Uc = 3 in both cases. Note 
that the Xcj decay constants may also depend on rUc- Obviously, the colour-octet 
contribution will also play an important role in the case of the Xhj decays. 

In the colour-octet contributions to the exclusive Xcj decays are estimated by 
calculating the hard scattering amplitude from the set of Feynman graphs shown in 
Fig.0 and convoluting it with the asymptotic pion wave function. The colour-octet 
and singlet contributions are to be added coherently. The Xcj —> tttt decay widths are 
given in terms of a single non-perturbative parameter k which approximately accounts 
for the soft physics in the colour-octet contribution. A £t to the data [^, ^ yields 
K = 0.16 GeV^ (with rUc = 1.5 GeV; Aqcd = 0.2 GeV) and the widths 


r(Xco( 2 ) ^ tt+tt ) = 49.85 (3.54) keV. 


(15) 


Gomparison with (^) reveals that the inclusion of the colour-octet mechanism brings 
predictions and data in generally good agreement. The value found for the parameter 
K has a reasonable interpretation in terms of charmonium properties and the mean 
transverse momentum of the quarks inside the pions. Results for the decays into 
pairs of uncharged pions are presented in Table 2. The quoted results for the pionic 
charmonium decays refer to a calculation within the standard HSA but similar good 
results are found when the modihed HSA is used pl|. The only soft parameter 

/Q\ 

appearing in the latter calculation is the octet-decay constant /} ^ of the charmonium 
state. 

Thus it seems that the colour-octet mechanism leads to a satisfactory explanation 
of the decay rates of the Xcj into two pions. Of course, that mechanism has to pass 
more tests in exclusive reactions before this issue can be considered as being settled. 
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Figure 2: Representatives of the various groups of colour-octet decay graphs. 



r(Xcj ^ 77*^77^^) [keV] 

BRiXcJ ^ 77 ^^ 77 ^) [%] 

i?2 — 0 

25.7 

1.81 

0.18 

0.091 

Exp. PDG H 

42± 18 

2.2 ±0.6 

0.31 ±0.06 

0.110 ±0.028 


Table 1: Decay widths and branching ratios of Xcj (colour-octet contributions 

included; rric = 1.5 GeV, Aqcd = 0.2 GeV). 


4. The 77 form factor 

Let us now turn to the case of the tt form factor and discuss the implications of the 
constraints on the pion’s wave function obtained from the 777 analysis. The leading 
twist result for the pion form factor can be brought into a form similar to (|^) 

= y |1 + + 0(7)]. (16) 


Ghoosing /iy = /i^ and using the value of determined in the leading twist 

analysis of the as well as a value of, say, 0.4 for 0*, one obtains the result 
0.097 GeV^/Q^ for to lowest order pQGD. That result is much smaller than the ad¬ 


mittedly poor experimental result [Q: Ft, = 0.35±0.10 GeV^/Q^. The as-corrections 


are too small to account for that discrepancy |^. The modihed HSA likewise pro- 


























vides too small a perturbative contribution p3. It is important to remember at this 


point that, formally, the perturbative contribution to the pion form factor represents 
the overlap of the large momentum tails of the initial and hnal state wave functions. 
But the form factor also gets a contribution from the overlap of the soft wave functions 
That contribution, frequently termed the Feynman contribution, is customarily 
assumed to be negligible already at momentum transfers as low as a few GeV^? Ex¬ 
amining the validity of that presumption by estimating the Feynman contribution 
from the asymptotic wave function (H), one hnds results of appropriate magnitude to 
hll in the gap between the perturbative contribution and the data of [Q. The results 
exhibit a broad flat maximum which, for momentum transfers between 3 and about 
15 GeV^, simulates the dimensional counting behaviour. For a wave function based 
on the Ghernyak-Zhitnitsky distribution amplitude, on the other hand, the Feynman 
contribution exceeds the data signihcantly i,@]- Large Feynman contributions have 
also been found by other authors O, ^ 


|35| . Thus, the small size of the perturbative 
contribution to the elastic form factor hnds a comforting although model-dependent 
explanation, a fact which has been pointed out by Isgur and Llewellyn Smith long 
time ago. Of course this line of reasoning is based on the assumption that the data 
of 


33] are essentially correct. 


A comment concerning the pion form factor in the time-like region is in order. In 
the standard HSA the predictions for the form factor in both the time-like and the 
space-like region, are identical. The experimental information on the time-like form 
factor comes from two sources, e’*'e“ —> 7r'''7r“ and J/T —>■ 7r’''7r“, which provides, 
to a very good approximation]! the form factor at s = Mj^^. Although the e"*" e“ 
annihilation data of Bollini et al. suffer from low statistics, they agree very well 
with the result obtained from the J/T decay. Gombining both the data, one hnds 
\Ft^\ = 0.93 ± 0.08 GeV^/s in the momentum transfer range between 2 and 10 GeV^, 
a value wich is roughly a factor of 3 larger than the space-like data. The modihed 
HSA can account for that large ratio of the time-like over space-like form factors 
although the perturbative contributions to the form factor in both the regions are too 
small. 

The structure function of the pion offers another possibility to test the wave func¬ 
tion against data. As has been pointed out in [Q the parton distribution functions 


are determined by the Fock state wave functions. Since each Fock state contributes 
through the modulus squared of its wave function integrated over transverse momenta 
up to Q and over all fractions x except those pertaining to the type of parton consid¬ 
ered, the contribution from the valence Fock state should not exceed the data of the 
valence quark structure function. As discussed in ||^, the asymptotic wave function 
respects this inequality while the Ghernyak-Zhitnitsky one again fails dramatically. 


^ At large the Feynman contribution is suppressed by l/Q^ as compared to the perturbative 
contribution. 


!! The contribution from cc transitions into the light quarks via three gluons cancel to zero if quark 
masses are neglected 371. 
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5. Summary 

The study of hard exclusive reactions is an interesting and challenging subject. 
The standard HSA, i.e. the valence Fock state contribution in collinear approximation 
to lowest order perturbative QCD, while asymptotically correct (at least for form 
factors), does not lead to a consistent description of the data. In many cases the 
predicted perturbative contribution to particular exclusive reactions are much smaller 
then the data. The observed spin effects do not find a comforting explanation. In 
some reactions agreement between prediction and experiment is found although at 
the expense of dominant contributions from the soft end-point regions rendering the 
perturbative analysis inconsistent. 

From a detailed analysis of the Try transition form factor it turned out that the pion 
distribution amplitude is close to the asymptotic form. Strongly end-point concen¬ 
trated distribution amplitudes are obsolete and the ostensible successes in describing 
the large momentum transfer behaviour of the pion form factor in the space-like re¬ 
gion and charmonium decays into two pions with such distribution amplitudes must 
therefore be dismissed. 

In view of these observations it seems that higher Fock state and/or higher twist 
contributions have to be included in the analysis of exclusive reactions. However, not 
much is known about them as yet. We are lacking systematic investigations of such 
contributions to exclusive reactions. The colour octet model for exclusive charmonium 
decays is discussed in this talk as an example of such contributions. 
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